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Abstract. Some of the means through which the possible presence of nearly 
. deconfined quarks in neutron stars can be detected by astrophysical observations of 

neutron stars from their birth to old age arc highlighted. 
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i 1. Introduction 
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Utilizing the asymptotic freedom of QCD, Collins and Perry [lj first noted that the dense 
cores of neutron stars may consist of deconfined quarks instead of hadrons. The crucial 
question is whether observations of neutron stars from their birth to death through 
neutrino, photon and gravity-wave emissions can unequivocably reveal the presence of 
nearly deconfined quarks instead of other possibilities such as only nucleons or other 
exotica such as strangeness-bearing hyperons or Bose (pion and kaon) condensates. 

2. Neutrino signals during the birth of a neutron star 
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The birth of a neutron star is heralded by the arrival of neutrinos on earth as confirmed 
by 1MB and Kamiokande neutrino detectors in the case of supernova SN 1987A. Nearly 
all of the gravitational binding energy (of order 300 bethes, where 1 bethe = 10 51 erg) 
released in the progenitor star's white dwarf-like core is carried off by neutrinos and 
antineutrinos of all flavors in roughly equal proportions. The remarkable fact that the 
weakly interacting neutrinos are trapped in matter prior to their release as a burst is due 
to their short mean free paths in matter, A ~ (cm) -1 ~ 10 cm, (here a ~ 10~ 40 cm 2 is the 
neutrino-matter cross section and n ~ 2 to 3 n s , where n s ~ 0.16 fm -3 is the reference 
nuclear equilibrium density), which is much less than the proto- neutron star radius, 
which exceeds 20 km. Should a core-collapse supernova occur in their lifetimes, current 
neutrino detectors, such as SK, SNO, LVD's, AMANDA, etc., offer a great opportunity 
for understanding a proto-neutron star's birth and propagation of neutrinos in dense 
matter insofar as they can detect tens of thousands of neutrinos in contrast to the tens 
of neutrinos detected by 1MB and Kamiokande. 

The appearence of quarks inside a neutron star leads to a decrease in the maximum 
mass that matter can support, implying metastability of the star. This would occur 
if the proto-neutron star's mass, which must be less than the maximum mass of the 
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Figure 1. Evolutions of the central baryon density ub, v concentration Y„, quark 
volume fraction \ an d temperature T for different baryon masses Mb- Solid lines 
show stable stars whereas dashed lines showing stars with larger masses are metastable. 
Diamonds indicate when quarks appear at the star's center, and asterisks denote when 
metastable stars become gravitationally unstable. Figure after Ref. [2] 



hot, lepton-rich matter is greater than the maximum mass of hot, lepton-poor matter. 
For matter with nucleons only, such a metastability is denied (see, e.g., [3]). Figured] 
shows the evolution of some thermodynamic quantities at the center of stars of various 
fixed baryonic masses. With the equation of state used (see [2] for details), stars with 
Mb < 1-1 M do not contain quarks and those with Mb ~ 1.7 M are metastable. 
The subsequent collapse to a black hole could be observed as a cessation in the neutrino 
signals well above the sensitivity limits of the current detectors (Figure [2]). 

3. Photon signals during the thermal evolution of a neutron star 

Mult i wavelength photon observations of neutron stars, the bread and butter affair of 
astronomy, has yielded estimates of the surface tempeartures and ages of several neutron 
stars (Fig. [3]). As neutron stars cool principally through neutrino emission from their 
cores, the possibility exists that the interior composition can be determined. The 
star continuosly emits photons, dominantly in x-rays, with an effective temperature 
T e ff that tracks the interior temperature but that is smaller by a factor ~ 100. The 
dominant neutrino cooling reactions are of a general type, known as Urea processes 
[1], in which thermally excited particles undergo beta and inverse-beta decays. Each 
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Figure 2. The evolution of the total neutrino luminosity for stars of indicated baryon 
masses. Shaded bands illustrate the limiting luninosities corresponding to a count 
rate of 0.2 Hz in all detectors assuming 50 kpc for 1MB and and Kamioka, 8.5 kpc 
for SNO, SuperK, and UNO. Shaded regions represent uncertanities in the average 
neutrino energy from the use of a diffusion scheme for neutrino transport in matter. 
Figure after Ref. [2]. 

reaction produces a neutrino or anti-neutrino, and thermal energy is thus continuously 
lost. Depending upon the proton-fraction of matter, which in turn depends on the 
nature of strong interactions at high density, direct Urea processes involving nucleons, 
hyperons or quarks lead to enhanced cooling compared to modified Urea processes in 
which an additional particle is required to conserve momentum. However, effects of 
superfluidity abates cooling as sufficient thermal energy is required to break paired 
fermions. In addition, the poorly known envelope composition also plays a role in the 
inferred surface temperature (Fig. [3]). The multitude of high density phases, cooling 
mechanisms, effects of superfluidity, and unknown envelope composition have thus far 
prevented definitive conclusions to be drawn (see, e.g., [5]). 

4. Mesured masses and their implications 

Several recent observations of neutron stars have direct bearing on the determination of 
the maximum mass. The most accurately measured masses are from timing observations 
of the radio binary pulsars. As shown in Fig. HI which is compilation of the measured 
neutron star masses as of November 2006, observations include pulsars orbiting another 
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Figure 3. Observational estimates of neutron star temperatures and ages together 
with theoretical cooling simulations for M = 1.4 M Q . Models and data are described in 
[6] . Orange error boxes (see online) indicate sources from which both X-ray and optical 
emissions have been observed. Simulations are for models with Fe or H envelopes, 
with and without the effects of superfluidity, and allowing or forbidding direct Urea 
processes. Models forbidding direct Urea processes are relatively independent of M 
and supcrfluid properties. Trajectories for models with enhanced cooling (direct Urea 
processes) and superfluidity lie within the yellow region, the exact location depending 
upon M as well as superfluid and Urea properties. Figure adapted from Ref. [7J. 

neutron star, a white dwarf or a main- sequence star. 

One significant development concerns mass determinations in binaries with white 
dwarf companions, which show a broader range of neutron star masses than binary 
neutron star pulsars. Perhaps a rather narrow set of evolutionary circumstances conspire 
to form double neutron star binaries, leading to a restricted range of neutron star masses 
[9]. This restriction is likely relaxed for other neutron star binaries. A few of the 
white dwarf binaries may contain neutron stars larger than the canonical 1.4 M value, 
including the intriguing case jTO] of PSR J0751+1807 in which the estimated mass with 
la error bars is 2.1 ± 0.2 M . In addition, to 95% confidence, one of the two pulsars 
Ter 5 I and J has a reported mass larger than 1.68 M 

Whereas the observed simple mean mass of neutron stars with white dwarf 
companions exceeds those with neutron star companions by 0.25 M , the weighted 
means of the two groups are virtually the same. The 2.1 M neutron star, PSR 
J0751 + 1807, is about 4a from the canonical value of 1.4 M . It is furthermore the 
case that the 2a errors of all but two systems extend into the range below 1.45 M & , so 
caution should be exercised before concluding that firm evidence of large neutron star 
masses exists. Continued observations, which will reduce the observational errors, are 
necessary to clarify this situation. 

Masses can also be estimated for another handful of binaries which contain an 
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Figure 4. Measured and estimated masses of neutron stars in binarie pulsars (gold, 
silver and blue regions online) and in x-ray accreting binaries (green). For each region, 
simple averages are shown as dotted lines; error weighted averages are shown as dashed 
lines. For labels and other details, consult Ref. [8]. 



accreting neutron star emitting x-rays. Some of these systems are characterized by 
relatively large masses, but the estimated errors are also large. The system of Vela X-1 
is noteworthy because its lower mass limit (1.6 to 1.7M Q ) is at least mildly constrained 
by geometry [12]. 

Raising the limit for the neutron star maximum mass could eliminate entire families 
of EOS's, especially those in which substantial softening begins around 2 to 3n s . This 
could be extremely significant, since exotica (hyperons, Bose condensates, or quarks) 
generally reduce the maximum mass appreciably. 

Ultimate energy density of observable cold baryonic matter 

Measurements of neutron star masses can set an upper limit to the maximum possible 
energy density in any compact object. It has been found [32] that no causal EOS has 
a central density, for a given mass, greater than that for the Tolman VII [TJ] analytic 
solution. This solution corresponds to a quadratic mass-energy density p dependence 
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Figure 5. Model predictions are compared with results from the Tolman IV and 
VII analytic solutions of general relativistic stucture equations. NR refers to non- 
relativistic potential models, R are field- theoretical models, and Exotica refers to NR 
or R models in which strong softening occurs, due to hyperons, a Bose condensate, or 
quark matter as well as self-bound strange quark matter. Constraints from a possible 
redshift measurement of z = 0.35 is also shown. The dashed lines for 1.44 and 2.2 M Q 
serve to guide the eye. Figure taken from Ref. [13]. 



on r, p = p c [l — (r/R) 2 ], where the central density is p c . For this solution, 

/M \ 2 

PcTvn = 2.5p cJnc - 1.5 x 10 16 (j^J gem" 3 . (1) 

A measured mass of 2.2 M would imply p max < 3.1 x 10 15 g cm -3 , or about 8n s . 

Figure [5] displays maximum masses and accompanying central densities for a wide 
wariety of neutron star EOS's, including models containing significant softening due 
to "exotica", such as strange quark matter. The upper limit to the density could be 
lowered if the causal constraint is not approached in practice. For example, at high 
densities in which quark asymptotic freedom is realized, the sound speed is limited to 
c/ \f?>. Using this as a strict limit at all densities, the Rhoades & Ruffini [15] mass 
limit is reduced by approximately l/\/3 and the compactness limit GM/Rc 2 = 1/2.94 
is reduced by a factor 3~ 1//4 to 1/3.8 I Til . In this extreme case, the maximum density 
would be reduced by a factor of 3 _1//4 from that of Eq. (pQ). A 2.2 M Q measured mass 
would imply a maximum density of about 4.2n s . 

5. Gravitational wave signals during mergers of binary stars 

Mergers of compact objects in binary systems, such as a pair of neutron stars (NS- 
NS), a neutron star and a black hole (NS-BH), or two black holes (BH-BH), are 
expected to be prominent sources of gravitational radiation [17]. The gravitational- 
wave signature of such systems is primarily determined by the chirp mass M c hi rp = 
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Figure 6. Physical and observational variables in mergers between low-mass black 
holes and neutron stars or self-bound quark stars. The total system mass is 6 M e and 
the initial mass ratio is q — 1/3 in both cases. The initial radii of the neutron star 
and quark star were assumed to be equal. The time scales have arbitrary zero points. 
Upper panel displays semi-major axis a (thick lines) and component mass M^s, Mqs 
(thin lines) evolution. Lower panel displays orbital frequency v (thick lines) and strain 
amplitude evolution. Solid curves refer to the neutron star simulation and dashed 
curves to the quark star simulations. Figure taken from Ref. |8J. 

(M l M 2 f /b {M 1 + Ma)" 1 / 5 , where M x and M 2 are the masses of the coalescing objects. 
The radiation of gravitational waves removes energy which causes the mutual orbits to 
decay. For example, the binary pulsar PSR B1913+16 has a merger timescale of about 
250 million years, and the pulsar binary PSR J0737-3039 has a merger timescale of 
about 85 million years [18J, so there is ample reason to expect that many such decaying 
compact binaries exist in the Galaxy. Besides emitting copious amounts of gravitational 
radiation, binary mergers have been proposed as a source of the r-process elements [19] 
and the origin of the shorter-duration gamma ray bursters [20] . 

Observations of gravity waves from merger events can simultaneously measure 
masses and radii of neutron stars, and could set firm limits on the neutron star maximum 
mass [2T1 122] . Binary mergers for the two cases of a black hole and a normal neutron 
star and a black hole and a self-bound strange quark matter star (Fig. [6]) illustrate 
the unique opportunity afforded by gravitational wave detectors due to begin operation 
over the next decade, including LIGO, VIRGO, GEO600, and TAMA. 

A careful analysis of the gravitational waveform during inspiral yields values for 
not only the chirp mass M c ^ rp , but for also the reduced mass MbhMns/M, so that 
both Mbh an d M^s can be found [23]. The onset of mass transfer can be determined 
by the peak in u, and the value of u there gives a. A general relativistic analysis of 
mass transfer conditions then allows the determination of the star's radius [22]. Thus 
a point on the mass-radius diagram can be estimated [21]. The combination h + uj^ 1 ^ 
depends only on a function of q, so the ratio of that combination and knowledge of 
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should allow determination of qj. From the Roche condition and knowledge of aj from 
ujf, another mass-radius combination can be found. 

The sharp contrast between the evolutions during stable mass transfer of a normal 
neutron star and a strange quark star should make these cases distinguishable. For 
strange quark matter stars, the differences in the height of the frequency peak and the 
plateau in the frequency values at later times are related to the differences in radii of 
the stars at these two epochs. It could be an indirect indicator of the maximum mass 
of the star: the closer is the star's mass before mass transfer to the maximum mass, the 
greater is the difference between these frequency values, because the radius change will 
be larger. Together with radius information, the value of the maximum mass remains the 
most important unknown that could reveal the true equation of state at high densities. 
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